The thermal Casimir effect for rough metallic plates 
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We propose a new theory of thermal Casimir effect, holding for the experimentally important 
case of metallic surfaces with a roughness having a spatial scale smaller than the skin depth. The 
theory is based on a simple phenomenological model for a rough conductor, that explicitly takes 
account of the fact that ohmic conduction in the immediate vicinity of the surface of a conductor is 
much impeded by surface roughness, if the amplitude of roughness is smaller than the skin depth. 
As a result of the new model, we find that surface roughness strongly influences the magnitude of 
the thermal correction to the Casimir force, independently of the plates separation. Our model, 
while consistent with recent accurate measurements of the Casimir force in the submicron range, 
leads to a new prediction for the not yet observed thermal correction to the Casimir force at large 
plates separation. Besides the thermal Casimir problem, our model is relevant for the correct 
theoretical interpretation of current experiments probing other proximity effects between conductors, 
like radiative heat transfer and quantum friction. 
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The Casimir effect [l[ provides physicists with the rare 
opportunity to investigate in a laboratory the physics 
of the quantum vacuum a subject that has recently 
attracted much interest for a number of reasons. One one 
hand, the possibility exists that the recently discovered 
accelerated expansion of Universe might be explained by 
a cosmological term originating from quantum vacuum 
effects [H . Another reason of interest comes from current 
experiments on non-newtonian forces at the sub-micron 
scale Q. As its is well known, Casimir forces are the 
dominant forces at this scale, and therefore determining 
them very accurately is indispensable in order to obtain 
experimental bounds on possible non-Newtonian forces 
Q. Finally we would like to mention that the Casimir 
effect, and more in general van der Waals forces M, are 
now finding new applications in nanotcchnology [7|. 

In the field of Casimir physics, an important and much 
controversial problem is that of determining how the 
Casimir force between two metallic bodies is modified 
by the temperature of the environment Apart from 
its intrinsic interest as a problem in theory of quan- 
tum electromagnetic (e.m.) fluctuations, addressing this 
problem is important because many experiments on non- 
newtonian forces at the submicron scale use metallic sur- 
faces at room temperature Q. We remark that, while 
the Casimir pressure has now been measured accurately 
[B|, [j| , the thermal contribution to the Casimir effect has 
not yet been detected as we write, and indeed there are 
at present several ongoing and planned experiments, to 
measure it 0. Surprisingly, the theory of the thermal 
Casimir effect for metallic bodies turned out to be rather 
complicated, and at present there exist several conflict- 
ing approaches that give widely different predictions for 
the magnitude of the effect 0. 



In this Letter we present a new theory of the thermal 
Casimir effect for conductors. The distinctive new fea- 
ture of the theory consists in the central role played in 
our approach by the surface roughness of the conductor. 
The presence of roughness is indeed a general feature 
of current Casimir experiments, which use metallic sur- 
faces presenting surface roughness with a characteristic 
amplitude ft, of a few nanometers P, H[. Until now, it 
has always been thought that for separations a between 
the surfaces much larger than the roughness amplitude 
h, surface roughness implies only a small correction of or- 
der {h/af 0]. We argue that this conclusion, while cor- 
rect for insulators, is not warranted for metallic surfaces 
at finite temperature, if the roughness h is smaller than 
the skin depth 5, because then account must be taken 
of the fact that surface irregularities impede ohmic con- 
duction in the immediate vicinity of the surface, a fact 
noted already by Pippard [l0( long ago. Using a sim- 
ple phenomenological model to account for this effect, 
we then find that independently of the separation a, the 
magnitude of the thermal correction to the Casimir force 
strongly depends on the surface roughness. While consis- 
tent with a recent precise measurement 5] of the Casimir 
pressure in the submicron range, our theory leads to a 
new prediction for the magnitude of the not yet mea- 
sured thermal correction to the Casimir force at large 
separations. We note that the application range of our 
model extends to other proximity effects originating from 
e.m. fluctutations, like radiative heat transfer ll| and 
quantum friction between closely space metallic surfaces 
[12j |. that are currently under intense investigation, but 
we shall not tough upon these issue in this Letter. 

All current approaches to the thermal Casimir effect 
are based on the theory of the Casimir effect developed 
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long ago by Lifshitz |13j , on the basis of Rytov's general 
theory of e.m. fluctuations 14J. Lifshitz-Rytov theory 



works very well in the case of media that can be char- 
acterized by electric and magnetic permittivity tensors, 
eij{to) and Hij{to) that depend only on the frequency, but 
are independent of the wave vector k, i.e. for media that 
exhibit only time-dispersion but not space-dispersion. 
According to this theory, the Casimir pressure P between 
two identical plane-parallel (possibly stratified) slabs at 
temperature T, separated by an empty gap of width a, 
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where the prime over the Z-sum means that the I = 
term has to taken with a weight one half, kj_ denotes the 
projection of the wave- vector onto the plane of the plates 
and qi = y/k\ + where & = 2i:k B Tl/h. The quan- 
tities r Q (i^;,kj_) denote the reflection coefficients of the 
slabs for a-polarization (for simplicity, we do not consider 
the possibility of a non-diagonal reflection matrix) , eval- 
uated at complex frequencies i^i . For the case of single- 
layer homogeneous slabs, with permittivity e(w), r a are 
the familiar Fresnel reflection coefficients. The above for- 
mula for the Casimir pressure is regularly used also in the 
case of metals, that can be characterized by a (complex) 
conductivity a(ui), in which case one has for the permit- 
tivity the expression e(uS) = l-\-&mo~{bS)/bJ. The difficul- 
ties found recently in the study of the thermal Casimir 
effect for metals originate from the fact that in order to 
evaluate the I = term of Eq. (JTJ one needs extrapolate 
the reflection coefficients r a to zero frequency. Different 
prescriptions for doing this have been proposed, leading 
to largely different predictions for the thermal correction 
to the pressure, depending on the resulting magnitude of 
the I = term. For example, the approach of Ref. 15 1 
adopts the Drude model 
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where Op is the plasma frequency and 7 the relaxation 
frequency accounting for ohmic conductivity. When this 
model is used, one finds that the I = term for TE 
polarization gives zero contribution, and the predicted 
thermal correction is much larger than what one obtains 
for the ideal case of perfect reflectors (corresponding to 
taking r 2 a = 1 in Eq. fl}). Recently, it has been shown 
that this model is inconsistent with the experiment [5j 
and moreover it has been noted that, in the case of per- 
fect lattices, the Drude model leads to thermodynamical 
inconsistencies at low temperature [l6| . On the contrary, 
the approach based on the plasma model 



e P (uj) = 1 - tt 2 P /L0 2 
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(augmented by a six-oscillator contribution accounting 
for interband transitions), has been shown to be consis- 



tent with the experiment [f|, and it is immune of ther- 
modynamic inconsistencies. 

All the above is for perfectly plane-parallel plates 
while, as we said earlier, real experiments involve rough 
surfaces, with amplitude h typically in the range of a 
few nanometers. In the existing literature, the effect of 
the slabs roughness has always been regarded as a geo- 
metric problem, involving as the only relevant parameter 
the geometric quantity h/a. For h <C a, the problem 
has been studied by perturbative means, for example the 
so-called Proximity Force Approximation (PEA) which 
corresponds to taking a suitable average of the Lifshitz 
formula result, across the surface of the plates {1]. The 
result of the PFA is that for h <C a, roughness gives only a 
small fractional correction of order (h/a) 2 to the Casimir 
force between two perfectly flat plates (see for example 
Ref. [U). Below we show that such a conclusion is not 
correct however, if the roughness profile varies rapidly 
along the surface of the conductor, over distances equal 
to the skin depth 8 of the e.m. fields, as it is usually the 
case in current Casimir experiments. 

Our analysis of the effect of roughness starts from the 
remarks made a long time ago by Pippard in his inves- 
tigations of the skin effect [lfj. Pippard noted that at 
frequencies in the GHz region, the presence of micro- 
scopic scratches on the surface of a metal can determine 
a significant increase of the surface resistance, by forcing 
the electric currents to flow along longer, distorted paths 
inside the metal. We considered that these remarks ac- 
quire even greater importance if the irregularities of the 
surface are smaller than the skin depth 8. If such a con- 
ductor is placed in a tangential alternating electric field, 
the induced ohmic currents in the immediate vicinity of 
the surface will be unable to follow the rapid distortions 
of the surface, with the result that barely any current will 
be found in the outer layer of thickness h of the conduc- 
tor. This is like water flowing inside a large pipe whose 
inner surface presents small bumps. The water contained 
in these bumps does not take part in the flow, and re- 
mains still. We were thus led to conceive that the surface 
layer of thickness h of a rough conductor does not behave 
like the deeper part of the conductor, and in fact, with 
regards to tangential electric fields, it rather resembles a 
dielectric, having a negligible real conductivity cr(u>). A 
detailed investigation of the response function of a rough 
metallic surface, according to the above picture, is likely 
to be a rather complicated problem, because when h is 
comparable to or less than the electron mean free path I 
at room temperature (for gold / is around 20 nm), spa- 
tial dispersion is expected to be relevant. In this Letter 
we shall content ourselves with proposing a simple phe- 
nomenological model, suggested by the previous qualita- 
tive considerations. We neglect the possible effect of spa- 
tial dispersion near the metal's surface, and we model a 
(thick) rough metallic plate as a plane-parallel two-layer 
system, consisting of a thick slab with permittivity e(w), 
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describing the inner part of the conductor, covered by 
a thin plane-parallel uniform layer of thickness h with a 
different permittivity e sur f(w), modelling its rough sur- 
face. For the bulk permittivity e(u>) we take the usual 
Drude model in Eq. @, while for the thin outer layer 
we take the plasma model: 



Esurf (w) = 1 - ^surf/^ 2 



with 
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where / is the fraction of the plane-parallel layer of width 
h really occupied by the metal. We now explain these 
choices. Modelling the rough surface as a uniform plane- 
parallel layer is reasonable, because in typical Casimir 
experiments at room temperature the wavelengths A ~ a 
and the skin depths S of the e.m. fields are both much 
larger than h, and therefore the e.m. fields effectively see 
the rough surface as flat. The choice of the plasma- model 
for e sur f appears as the simplest one, because the plasma 
model is what one obtains from the Drude model, after 
one sets to zero the real part of the conductivity a, in 
such a way that the surface layer behaves as an insulator, 
as we argued it is the case. The chosen value of Sl 2 urf in 
Eq. {5J depends on the fact that the average density 
"surf of free electrons within the rough layer is smaller 
by a factor / than the bulk value n. Upon recalling 
that the square of the plasma frequency is related to the 
density n of free electrons by the well known relation 
£1 2 = Aime 1 /m, where e and m denote the electron and 
(effective) mass, respectively, one arrives at Eq. (0. In 
principle, for a given rough surface, it should be possible 
to compute the values of h and /, starting from a detailed 
microscopic theory. In practice, we shall regard h and 
/ as phenomenological parameters to be determined by 
comparison with experimental data. 

The two layer structure for a rough metallic plate is 
what distinguishes our model from the models that have 
appeared in the recent literature, which all use a single 
dielectric function e(to) for the entire plate, the various 
approaches differing in the choice of e(ui). It is clear that 
for sufficiently large h our model is expected to reproduce 
the results of the plasma model used in Ref. [5j , while for 
h — ► it should resemble the Drude model of Ref. [la]. 
To see exactly how things go, we now compute the ther- 
mal Casimir force between two identical rough plates, 
using our model. This is done simply by substituting 
into Lifshitz formula Eq. ^ the appropriate reflection 
coefficients r a for our two-layer model of a rough plate: 
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fcj, and are the usual Fresnel re- 
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FIG. 1: Experimental data for the reduction factor 77 from 
Ref. [f| (diamonds). The inset shows a magnified view of 
the picture for small separations. The solid and dotted lines 
are fits using our theory, with h = 11 nm and / = 0.9, for 
T = 300 K (solid line) and T = K (dotted line). The dot- 
dashed line is for the Drude model approach at T = 300 K. 
Distances d are in microns. See text for explanations 



we shall plot the Casimir pressure as a function of the 
average separation d between the plates, that is related 
to a as d = a + 2h(l — /). Indeed d is a more meaningful 
quantity to consider, because it coincides with the plates 
separation that is usually measured in the experiments. 

As a check of our theory, we tried to fit a recent (indi- 
rect) accurate measurement of the Casimir force between 
two plane-parallel gold surfaces, using a micromachincd 
oscillator, reported in Ref. The experimental data 
are shown in Fig. 1, in terms of the reduction factor 
jj = P/Pid, where P id = (7r 2 ^c)/(240d 4 ) is the Casimir 
pressure for two ideal metallic plates, at zero tempera- 
ture. In our computations, we took hQp = 8.9 eV and 
h"f — .0357 eV, which are the values used in Q. As the 
most accurate measurements were performed at rather 
small separations, around 200 nm, interband transitions 
of core electrons give a sizable contribution to the per- 
mittivity, and we therefore added both to e and e sur f the 
six-oscillator expression cib(w) accounting for interband 
transitions, that was used in Ref.0]. Fig.l shows two fits 
with our model for h = 11 nm and / = 0.9, for T = 300 
K (solid-line) and T = (dotted line). We see that both 
curves are consistent with all the measurements in the 
entire range of distances, ranging from 162 nm to 746 
nm, with errors well inside the 95 confidence interval re- 
ported in Ref. [||. The values for h and / obtained from 
the fit are nicely consistent with the degree of roughness 
of the Au surfaces reported in Ref. || , which appears as 
encouraging in favor of our model. As we see, the preci- 
sion of the measurements is not sufficient to discriminate 
the T = and T = 300 K curves, and therefore this 
experiment does not detect the thermal contribution to 
the Casimir force. We remark that the data rule out 
the Drude-approach of Ref. 



15j (dot-dashed line), while 



flection coefficients for the interface ij. In what follows they are consistent with the plasma- model approach ad- 
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FIG. 2: Plots of the reduction factor r\ for Au plates at T = 
300 K, as a function of the average separation d (in microns). 
The solid and dashed curves are for our theory, with h = 11 
nm and / = 0.9 (solid curve), h — 2 nm and / = 0.5 (dashed 
curve). The dotted curve is for the plasma- model of Ref. [a] 
and the point-dashed curve for the Drude model of Ref. (la ]. 



the scale of the skin depth, the usual approximations 
based on the PFA are expected to be correct, and then 
roughness should only give a small correction of order 
(h/a) 2 . We do not know yet if our two-layer model is free 
from the thermodynamical inconsistencies at low temper- 
ature that plague the Drude model [l6j]> an d we plan to 
investigate this important problem in a future work. An- 
other important remark is that the model presented in 
this Letter can be used to obtain predictions for other 
proximity effects, originating from thermal e.m. fluctua- 
tions, like for example the power S radiative heat trans- 
fer between two metallic plates at different temperatures, 
separated by an empty gap. Investigating this problem 
would be very interesting, because that would provide 
us with another check of the correctness of our model 
for a rough metallic surface. This subject also will be 
addressed in future work. 



vocated in Ref. [5[, and in fact for this experiment our 
model gives results that are very close to those obtained 
by means of the plasma model. 

It is now interesting to consider how the magnitude 
of the thermal correction depends on the roughness at 
larger separations, and in Fig. 2 we show plots of the 
reduction factor for different degrees of roughness. The 
curves were computed for gold, with Mlp = 8.9 eV, and 
= .0357 eV, for T = 300 K. The dotted and dot- 
dashed line correspond, respectively, to the (generalized) 
plasma approach of Ref. [f| and to the Drude approach of 
Ref.[15j. The solid and dashed lines were computed us- 
ing our model, with / = 11 nm, / = 0.9 (solid) and I = 2 
nm, / = .5 (dashed). Fig. 1 lends itself to several im- 
portant comments. First of all we see that, according to 
our model, surface roughness does influence strongly the 
Casimir pressure also at large separations, contrary to the 
PFA picture. Second, neither the plasma model, nor the 
Drude model are generally valid. The plasma model may 
be expected to be valid only for "very" rough surfaces, 
having a roughness amplitude h around ten nm, like those 
used in the experiment in Ref. [5j. On the contrary, the 
Drude model should be valid for very smooth surfaces, 
with a roughness amplitude h well below two nm. For 
intermediate values of the roughness amplitude, neither 
theory should be adequate and perhaps the simple model 
proposed in this Letter provides a better approximation. 

In conclusion, we have presented a new phenomeno- 
logical theory for the thermal Casimir effect between two 
rough metallic plates. According to this theory, surface 
roughness influence strongly the thermal correction to 
the Casimir force at all separations a. The theory is 
expected to be valid for surfaces with a roughness hav- 
ing a characteristic spatial scale h below the skin depth 
<5, which represents the standard experimental situation. 
For smoother surfaces, having slow profile variations on 



* Bimonte@na.infn.it % 

[1] M. Bordag, U. Mohideen and V.M. Mostepanenko, Phys. 
Rep. 353, 1 (2001). 

[2] P. Milonni, The Quantum Vacuum (San Diego, CA: Aca- 
demic, 1994). 

[3] P.J.E. Peebles and B. Rata, Rev. Mod. Phys. 75, 559 
(2003). 

[4] E. Fischbach and C.L. Talmadge The Search for Non- 
Newtonian Gravity (New York: AlP/Springer, 1999). 

[5] R.S. Decca et al, Eur. Phys. J C 51, 963 (2007). 

[6] V. A. Parsegian, Van der Waals Forces (Cambridge Uni- 
versity Press, Cambridge, England, 2005). 

[7] Capasso F, Munday J N, Iannuzzi D and Chan H B 2007 
IEEE J. Sel. Top. Quant. Electron. 13 400 

[8] S.K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997); U. Mo- 
hideen and A. Roy, ibid. 81, 4549 (1998); G. Bressi, 
G. Carugno, R. Onofrio and G. Ruoso, ibid. 88 041804 

(2002) ; R.S. Decca, D. Lopez, E. Fischbach, and D.E. 
Krause, ibid. 91, 050402 (2003). 

[9] M. Brown-Hayes et al., J. Phys. A 39, 6195 (2006). 

[10] A. B. Pippard, Proc. Roy. Soc. A 191, 385 (1947) 

[11] D. Polder and M. Van Hove, Phys. Rev. B 4, 3303 
(1971); A.I. Volokitin and B.N.J. Persson, Phys. Rev. 
B 69, 045417 (2004); A. Kittel et al., Phys. Rev. Lett. 
95, 224301 (2005). 

[12] J.B. Pendry, J. Phys.: Condens. Matter 9, 10301 (1997). 

[13] E.M. Lifshitz, Sov. Phys. JETP 2, 73 (1956); E.M. Lif- 
shitz and L.P. Pitaevskii, Landau and Lifshitz Course 
of Theoretical Physics: Statistical Physics Part II 
(Butterworth-Heinemann, 1980). 

[14] S.M. Rytov, Theory of Electrical Fluctuations and Ther- 
mal Radiation, Publyshing House, Academy os Sciences, 
USSR (1953). 

[15] M. Bostrom and B.E. Sernelius, Phys. Rev. Lett. 84, 
4757 (2000); B.E. Sernelius, ibid. 87, 139102 (2001). 

[16] V.B. Bezerra, GL. Klimchitskaya and C. Romero, Phys. 
Rev. A 65 012111 (2001); B. Geyer, GL. Klimchitskaya 
and V.M. Mostepanenko, ibid. A 67 062102 (2003); V.B. 
Svetovoy and M.V. Lokhanin, Phys. Rev. A 67, 022113 

(2003) . 



